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Abstract 

The rheological behavior of metallic alloys containing both solid and liquid 
phases is investigated in the low solid fraction range (<50%). This behavior de- 
pends on both the solid fraction and the shear rate. The concept of Effective Vol- 
ume Fraction (EVF) is used to decorrelate the influence of these two parameters. 
At high shear rate the slurry behaves like a suspension of hard spheres, whereas at 
lower shear rate, particles tend to aggregate in clusters, entrapping liquid and thus, 
increasing the EVF and the viscosity. A lattice model is introduced to simulate 
the aggregation / break-up processes within a slurry under shear. When the steady 
state is reached, the entrapped liquid fraction is calculated, leading to a viscosity 
estimation. Simulation results for the viscosity and 3D cluster structure are in good 
agreement with experimental results. 

1 Introduction 

Semi-solid slurries are characterized by the coexistence of solid and liquid phases. They 
are usually observed in alloys with two or more constituents for temperatures between 
the solidus and liquidus lines. The rheological properties of this dual phase state is of 
interest both for casting and for metal forming operations known as thixoforming 
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When submitted to shear, the steady state viscosity decreases as the applied shear rate 
increases, reaching values of order 100 mPa.s @. 

This behavior is called rheofluidization and it is usually explained by the interaction 
between the solid particles. It is a very general phenomenon observed in suspensions, 
metalic alloys in the semi-solid state, colloids, latexes, etc. Their characteristic behavior 
is qualitatively interpreted as a competition between the aggregating (Coulomb attraction, 
Wan Der Waals forces, surface forces...) and break-up forces (shear). At low shear rate 
(typically lower than 1 s" 1 ), individual particles can aggregate in "clusters", being able 
to form a more or less rigid network and the slurry is considered as a solid. At very high 
shear rate, the motion of the particles prevents the particle/particle bonding and leads to 
a more dispersed suspension behaving more like a fluid. 

The concept of effective volume fraction introduced by Quemada || allows to relate 
the cluster characteristics to the viscosity of the mixture. The aim of the present paper 
is to build on this idea by considering the dynamics of cluster formation, and to derive 
a phenomenology for non-newtonian behaviour of semi-solid slurry, which is based on an 
understanding of the elementary phenomena governing the cluster dynamics. 

The rheofluidization has been extensively studied (see for instance the work of Que- 
mada @]). Figure [l| shows a typical shear rate dependence of the viscosity for an in- 
teracting particles suspension. Shear thickening that may occur for higher values of the 
shear rate || is not considered. The behavior at high shear rate when particles are well 
separated has been accurately described by Krieger H: 



-2.5$, 



$ being the solid fraction of the suspension and <&m its maximum value (close packing). 
Below a given value j s , an increase of the viscosity is observed. This can be interpreted 
as an increase in the effective volume fraction $ e // ; which takes into account both the 
solid and the entrapped liquid. $ e // replaces $ in equation [j]. 

Jeffrey and Acrivos 0] underlined the notion of structure resulting from aggregation 
between solid particles: for the same solid fraction $, different structure can lead to 
different viscosities. The structure depends itself on the shear rate. This statement is 
easily understandable for the two extreme limits of the shear rate: "percolating network" 
for low shear rate and "disperse suspension" for high shear rate. Between these two 
domains, the structure of diphasic solutions is not fully understood due to several reasons: 

• Mean-field modeling fails to describe the dynamics of such slurries because it does 
not take into account the interaction between solid particles. 

• 2D observations are useful but not sufficient because they do not show the real 
state of aggregation. In particular, 3D connectivity of the solid particles cannot be 
easily revealed. Serial cutting || allows 3D reconstruction, but is only tractable for 
relatively small zones. 

• Difficulty of 3D observations: Small Angle Scattering is efficiently used with 
colloids, but gives only a characteristic length-scale of the suspension. Moreover, it 
is limited to small solid particles (less than a few microns). A new observational 
technique is under development: synchrotron-radiation microtomography allowing 
3D investigation by means of phase contrast of a relatively large part of a sample 



with a resolution of a few micrometers [T(| 



In a previous work |Tl| , an analytical estimation of the characteristic radius of clus- 
ters was derived considering that aggregation and coalescence lead to spherical compact 
clusters. This radius was found to decrease with the shear rate as 7 _4//7 . In this pa- 
per we discuss the aggregation and break-up phenomena leading to more or less opened 
structures neglecting, this time, the densification of the resulting clusters. 

Different approaches has been used to simulate the rheological behaviour of colloidal 
suspensions. Dynamic simulation tools such as nonequilibrium Brownian dynamics [T2[ 
and Stokesian dynamics |13|] have provided insights into the influence of the various col- 
loidal forces on the microstructure and the effect of this microstructure on suspension 
rheology. Non-dimensional structural parameter has been introduced to predict the rhe- 
ology of aggregated sediment suspension |14| and semi-solid slurries [[L5| . 



Based on simple physical assumptions, we propose in this paper an original computer 
simulation of the microstructure of a semi-solid metallic alloy submitted to shear. First, 
we discuss the parameters governing the structure, then we use a simple approach to 
relate the microstructure with rheological properties. 



2 Structure and viscosity of sheared suspension 

The role of the solid fraction, $, in the rheology of sheared suspensions has been exten- 



sively studied (see for instance the review of Rutgers |16|). The most used law for the $ 
dependence of the viscosity is the Krieger phenomenological law. It has the advantage of 
pointing out the concept of maximum packing fraction (f>M- 

The role of the shear rate is less clear. We know that it tends to break-up the clusters: 
depending on the applied shear rate and the solid fraction. The possible structures can 
be schematically grouped in four different classes (see figure 0): 

• At low solid fraction ($ < $m), the structure depends on the shear rate. At low 
shear rate (1), very little break-up occurs, leading to a 3D interconnected network. 
At high shear rate (2), the break-up process dominates leaving the system as a 
dispersed suspension (individual particles are well separated). In between those two 
domains (3), the aggregation and break-up processes counterbalance, leading to a 
suspension of clusters with 7 dependent characteristic size and shape. 

• At high solid fraction ($ > $m), close-packing is observed for any applied shear 
rate. The shear is then localized in particular planes. 

The frontiers delimiting the different domains in figure |^ should be seen more as 
transition zones than strict boundaries. 

In the dispersed suspension region the viscosity can be accurately described with a 
model developed for suspensions of hard spheres (see equation |l]). In the 3D intercon- 
nected network region the deformation mechanisms are totally different involving much 
more the solid than the liquid. Pseudo plastic models leading to a power law, rj = raj 11-1 , 
describes successfully the rheology of the two-phased material (see reference 1 17]) . 

In the intermediate region, Quemada || introduced the concept of Effective Volume 
Fraction (EVF, $ e // in equation ||) to replace the solid fraction $ of equation [j], turning 
it into: 

V = VoU-^) (2) 



M 



The EVF takes into account the entrapped liquid not involved in the hydrodynamic 
flow. $ e // is then the sum of the real solid fraction $ and the part of the liquid which does 
no longer take part in the liquid flow. It could be, for example the entrapped liquid in the 
middle of a solid aggregate. Note that liquid does not need to be completely surrounded 
by solid to be entrapped (see section f|). On this basis the structure of the suspension 
is re-injected into equation p] through its correlation with the EVF. For a suspension of 
compact clusters, the EVF would be equal to the solid fraction, while for a more open 
structure (containing more entrapped liquid) it could be much higher than $ leading to 
higher values of viscosity. The maximum packing fraction &m characterizes a geometrical 
compactness. It is a function of the size and shape distribution of the individual solid 
particles and does not depend on 7. The EVF takes into account the dependence of the 
microstructure as a function of the applied shear rate. The structure could have been 



introduced in a shear dependent packing fraction $m HU, instead of $. Both approaches 
lead to the same results. 

Ito and Flemings || used this concept to interpret experimental results on semi-solid 
slurries. After a time consuming micrography analysis on many planes, they managed to 
give a 3D reconstruction of clusters resulting from stirring at steady state. They gave 
an estimation of the entrapped liquid fraction and plotted the measured viscosity versus 
EVF for slurries solidified under different shear rates. All the experimental points were 
scaled back on a unique master curve. 

The critical point of this analysis is the estimation of the structure dependent EVF, as 
a function of the shear rate. In the following section, we propose an approach predicting 
the structure evolution of a semi-solid slurry submitted to shear. 



3 The steady state of sheared suspension: competi- 
tion between aggregation and break-up 

When a semi-solid slurry is left at rest, two different mechanisms occurs: 



Oswald ripening |19] and spheroidization tend to minimize the solid/liquid surface 



energy by narrowing the size distribution of aggregates, and smoothing the solid par- 
ticles surfaces. These phenomena are diffusion limited with a characteristic length 
of the order of the individual particles radius a^. A characteristic diffusion time 
could be approximated by: a^/D m 2 s, where D « 5 x 10 -9 m 2 s~ 1 is the diffusion 
coefficient of the solute in the liquid and a® ~ 100 ^m. 

Particles aggregation occurs to lower the liquid/solid surface energy. The driving 
force of aggregation is the difference between twice the solid/liquid surface energy 
and the solid/solid surface energy: Act = Iosl — &b- The aggregation kinetic is 



limited by the collision frequency which is estimated in [^0] for a sheared suspension 
of hard spheres as: f c = {8/n) x $7. In the intermediate region of figure |2], the 
characteristic time between two collisions (l// c ) ranges from 1 s to 1 ms when 7 
ranges from 10 to 1000 s -1 . 

In the present paper, we are interested in a shear rate domain where the kinetic of 
aggregation is much more rapid than the kinetic of Oswald ripening (see above). We will 
then neglect Oswald ripening. Indeed, in a relevant time for diffusion limited mechanisms, 
aggregation induces more drastic changes in the microstructure than Oswald ripening. 



The aggregation mechanism is collision induced. We are interested in the physical 
mechanisms taking place during the contact between two particles. Two spherical objects 
of radius ao aiming at one another in a shear flow with shear rate 7 will meet during a 



mean contact time t c approximated by Adler [21| for the limit of small solid fraction: 



r c = \r x (3) 

During the contact, if a favorable crystallographic orientation is encountered, a rigid 
neck will be built between the two particles (see figure |3|). We now use a result derived in 
reference |15|] concerning the growth rate of a neck between two solid spherical particles 
of radius ao (see figure |3[). If x is the neck radius and Ty the neck formation time after 
the particles collision, Tf follows: 

where A is a constant dependent on thermo-physical properties of the material (see 



reference [11] for more details). Assuming the neck builds up during the contact time r c 



we have r c = r/, leading to an estimation of the neck size: 



x = 77T 7 5 «o ( 5 ) 



\2A, 

The shear will tend to break the aggregates during collision. As collision of two clusters 
of size AL occurs their mean relative velocity is AL/2 x 7. Thus, the kinetic energy E c 
available for rupture is: 

1 / A T \ 2 

p, V, 7 and AL are the mass density, the volume, the shear rate and the size of the 
cluster in the normal direction of the rupture plane, respectively. E c has to be compared 
with the energetic cost of the plastic rupture (classically considered to be 1000 times 



larger than the fragile rupture energy |f22f ): 



E r = n x irx x lOOOAa (7) 

n, and Aa being the number of broken necks of size x (defined in equation [5]) in the 
fracture plane, and the energy cost due to liquid / solid surface creation, respectively. 

If 7 is small, the contact time r c is long and necks will be large. Moreover, kinetic 
energy will be low preventing from any break-up of the resulting cluster. This situation 
will lead to the formation of a rigid 3D network. For large 7 the contact time is short 
and the kinetic energy is large enough to prevent any aggregation. This leads to a dis- 
persed suspension. As both aggregation and break-up are collision induced, the collision 
frequency influences only the kinetic of aggregation / break-up, but it does not influence 
the equilibrium between these two mechanisms. 

In the present study, we deal with " liquid driven deformation" . This means that the 
shear rate is sufficiently high to allow the competition between aggregation and break-up, 
the deformation being allowed by the free liquid between the clusters. 

From these simple considerations the probabilities of aggregation and break-up be- 
tween particles and clusters of particles will be derived and used to simulate cluster 
dynamics. 



4 Lattice model simulation of cluster structure 



In a real system containing « 10 5 particles, one would have to treat « 10 3 clusters in 
order to be representative of the distribution. This is untractable with standard computer 
power. We are interested here in steady state interaction, which means a dynamical equi- 
librium between cluster aggregation and cluster fracture. We assume a sort of "ergodic" 
hypothesis, namely that in steady state, the size distribution of a population of clusters 
at a given time, is equivalent with the size distribution of a given cluster along time. Our 
main simplifying approach based on the above hypothesis is in considering only one rep- 
resentative cluster that can either aggregate with its copy, or can be broken in a random 
plane. The time evolution of this cluster will approximate the ensemble properties of the 
equilibrium semi-solid slurry. 

The cluster is made up of a cubic arrangement of connected spheres (elementary 
particles of radius a ). It is stored in a 3D matrix (cubic lattice), where stands for free 
liquid, 1 stands for solid and 2 stands for entrapped liquid. The size of the matrix unit 
cell (2oq) is chosen so that the elementary particle is inscribed into this cubic unit cell. 

Aggregation 

The aggregation probability P a per collision is assumed to be the probability q e for two 



particles to encounter with a favorable crystallographic orientation. As in reference [15 
we consider that only low angle grain boundary (< 0.25 rad) give rise to a non wetted 
grain boundary, leading to: q e = 0.02. This condition is equivalent to: 



P a = q e 



Break-up 



The rupture probability per collision P r should be a monotonous increasing function of 
the ratio of the available kinetic energy over the energetic cost of the rupture (E c /E r ). 
At zero kinetic energy it should be zero, and for E c = E r it should converge to 1. The 
most simple form satisfying the above criterion is the simple P r = E c j E r approximation. 
A fracture plane is randomly chosen. The differential kinetic energy introduced in [ITT 



is 



given by equation |[ The volume V of the cluster is taken as (e + s)(2a ) 3 with e and 
s being the number of entrapped liquid voxels and the number of solid particles in the 
cluster, respectively. The energetic cost of the rupture is given by equation 0. We then 
have the rupture probability per collision: 

21 2 , . / at \2 , 12 

_E C _ 2 paj /25\" + ^ 5 (9) 

r E r i25ttAct \2AJ n 1 ' 

The first group of terms in ^| depends on material properties (p, Aa, A) and the initial 
state (do), whereas the second member is calculated at each step of the simulation. 

During cluster /cluster collisions, aggregation or rupture may occur. However, most of 
the time, neither aggregation nor break-up occur and the two clusters are left as before. 
As we deal with the steady state there is no need to consider these collisions with no 
effect. The only parameter is then the ratio P a /P r . 
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Through equation [TO], X incorporates the material properties and the initial state. For 
an Al-6.5wt%Si alloy with a 50 jum globular structure, K = 5 x 10 8 s -12 / 5 (see Table p. 

For each step, the one-cluster algorithm randomly selects a rupture plane and cal- 
culates the break-up probability (aggregation probability is fixed to q e ). A number is 
randomly generated between and P a + P r : 

• If it ranges between and P r fracture is processed by deleting all the solid particles 
belonging to the rupture plane. One of the resulting clusters is randomly selected 
and stored in the matrix M. 

• Else, the cluster is duplicated in a matrix N, randomly rotated, and the two matri- 
ces, M and N, are then aggregated allowing interpenetration. The resulting cluster 
is stored in the matrix M. 

Aggregation and rupture mechanisms are illustrated in figure [| and a simplified algo- 
rithm of the numerical simulation is shown in figure [|. 

The number of solid particles s, the number of entrapped liquid voxels e and the 
gyration radius R are extracted from the matrix M, leading to the calculation of the 
effective volume fraction. The gyration radius is defined as the mean distance between 
the cluster center of mass and all its solid particles. Averages (R, s, e) are calculated over 
the total step number. Since the mean "life-time" between two collisions is assumed to 
be independent of the cluster radius, there is no need to weight the averages. 

One of the crucial point of this simulation is the evaluation of the entrapped liquid 
in the cluster. Indeed, this effect governs the effective volume fraction $ e // that will be 
introduced in equation to predict the viscosity. It is calculated at each step and defined 
as follows: a liquid voxel is entrapped if 4 or more of its 6 principal directions (cubic lattice) 
(Ox, —Ox, Oy, —Oy, Oz, —Oz) hits any solid particle of the cluster. If only 3 directions 
out of the 6 point to a solid particle, the liquid is entrapped if the averaged distance 
between the considered liquid voxel and the 3 intersected solid particles is smaller than 
the cluster gyration radius R. In other terms, if d\, d 2 and d 3 are the distances between 
the liquid voxel and the nearest solid particle in the 3 considered directions, the liquid is 
entrapped if d\ + d\ + d\ < R? . For a better understanding, this procedure is visualized 
in figure for the 2D case. 



5 Numerical validation: 

The simulation starts with a cluster constituted of one solid particle. It usually growths 
until it reaches a large enough size to be broken. A typical time evolution of the cluster 
gyration radius R is shown in figure [?| for a shear rate of 500 s" 1 . Although the gyration 
radius R fluctuates strongly in time, the distribution of R is stable. The numerical 
convergence is reached when fluctuations of the distribution are less than 10%. Figure § 
displays the cluster gyration radius distribution after various computation steps. Steady 
state is reached after « 5000 steps. 



6 Results 



The result of the simulation will be presented as follows: 

1. 3D cluster shape at different shear rates will be depicted. 

2. The size distribution of clusters for different shear rates will be plotted and compared 
with an analytical model available in the literature. 

3. The 7 dependence of the mean radius R of the clusters will be discussed. 

4. The parameter $ e yj/<3> will be computed as a function of the shear rate giving a 
more precise structure map than the one given in figure |[ 

5. Finally, the viscosity rj will be given as a function of the shear rate and the solid 
fraction. 



Cluster shape 

Figure |9] shows a cluster of radius R with its entrapped liquid for two different shear rates. 
For small j, aggregation is dominant leading to a more open structure entrapping more 
liquid. 



Cluster size distribution 

Figure ^ exhibits normalized radius distributions for different shear rates. Size distribu- 
tion of particles can be accurately described by a "log-normal" law, frequently encountered 
in fragmentation problems (see reference ||23|1 ). Since these experiments are ensemble dis- 
tributions, whereas our simulation deals with time evolution, this similarity gives some 
confidence into the used "ergodic"-type hypothesis. 

The normalization brings the simulated distribution for various shear rates to a single 
master curve. The size dispersion is then proportional to the mean size of the clusters: 
lower is the shear rate, larger is the mean cluster radius and larger is the size dispersion. 
For low shear rate, the cluster can explore either small and large sizes, whereas, for high 
shear rate, only small sizes are possible. 

In figure [10], one can compare the simulated distribution with the mathematical model 



of Takajo [^4| often used as a model for aggregation processes. This model based on 
homogeneous coalescence frequency gives a normalized size distribution at steady state. 
These two different approaches lead to comparable profiles, giving some confidence into 
the general basis of our simulation. 



Mean cluster radius vs. shear rate 

The obtained mean gyration radius values are plotted as a function of the shear rate on 



figure [LI] The mean gyration radius R decrease with increasing shear rate following a 
power law R oc 7~ ar . This dependence is comparable with the power law estimated for 
the compact cluster scenario |TT| : R oc A f~ os . It can be observed that, above a given 7 
only elementary particles remain, hence there is no more evolution of R. 

As the physical parameter K depends strongly on the size of the individual particles 
a Q (see equation [TtJ) , computations are presented for two different reasonable values of a . 
For the smaller elementary particle size, cluster are larger at some given shear rate and 
solid fraction. 



Effective volume fraction vs. shear rate 



Figure [12] displays the calculated ratio $ e ///$ as a function of 7 for two different elemen- 
tary particle sizes. For high 7, (high enough to break all the necks), there is no entrapped 
liquid ($ e ^//$ = 1). For small 7, more connected structures leads to high fraction of 
entrapped liquid ($ e ///$ « 2). 

In figure we can identify the three different domains introduced in section |2| (fig- 
ure |2|) . 



For high 7 values (> 7,), the dispersed suspension-like behavior is observed. j s is 
reached when the mean radius of the clusters is equal to the radius of the elementary 
solid particles, that is to say when the effective volume fraction is equal to Note 



that an approximation of % can be obtained equaling P a and P r of equation 10. We 



find for Al-6.5wt%Si with initial particles size ao = 100 fim: 7 S = 1300 s . This 



value is clearly in accordance with the transition domain of figure |T2 



For low 7 values (< j n ) a rigid network sets in. The estimation of 7„ depends on 
the gelation mechanism: 

— experimentally, the apparent rigid network will appear when the mean radius 
of the clusters reaches the characteristic size of the measurement apparatus 
(Couette Rheometer for example). 

— in theory, the apparent rigid network will appear when EVF = $m : the whole 
volume of the sample is invaded by a unique cluster. 

7„ is found to be $ dependent. We can observe that the simulation leads to a more 
realistic description of the structure (figure [13]) than in figure |2|. Indeed, lower is the 



solid fraction, higher is the liquid fraction involved in the hydrodynamic flow for 
some given 7 and aggregation state characterized by e and s. 

• In the intermediate range we have a suspension of clusters with decreasing size and 
entrapped liquid fraction as 7 increase. 

Viscosity vs. shear rate 

For the calculation of the viscosity, the entrapped liquid is taken into account trough (a) 
the mean number of entrapped liquid voxels e, and (b) the amount of entrapped liquid 
within a cubic voxel containing a solid sphere, which depends on the mean number n n of 
solid neighbors of a particle over the 26 possible neighbors in the cubic lattice. The liquid 
fraction within a unit cell is 0.42 (difference of volume between the cube and its inscribed 
sphere). The viscosity follows: 



V = Vo (l 



*.// = * (1**1'+^ 



(ii) 



where e and s are the mean number of entrapped liquid sites and the mean number of 
solid sites in the cluster, respectively. With a suspension of hard spheres, $m would have 
been 0.65, but as we deal with polydisperse and deformable solid particles, it is chosen to 
be 1. Indeed, the liquid phase has been proven to remain connected for solid fraction up 
to 0.8. 



These equations lead to figure [b| where viscosity is plotted versus the solid fraction. 
The viscosity exhibits the typical profile of a semi-solid mixture @. The viscosity in- 
creases with the solid fraction until the solid fraction reaches a threshold value where the 
suspension behaves more like a solid than a liquid. Higher is the shear rate, higher is this 
threshold value. 



7 Comparison with real microstructure and experi- 
mental viscosity measurement 

7.1 3D microstructure visualization 

The advantage of such a simulation is the ability of viewing 3D structures. It is of interest 
to compare them with real structures resulting from serial cutting and micrography anal- 



ysis made by Ito ||. Figure [15] shows experimental and simulated 3D cluster structures. 
The physical parameter K (see equation |10|) has been calculated for the Al-6.5wt %Si alloy 
with a elementary particle radius ao — 50/im. Note that the cluster exhibited by Ito is 
not representative of the population. Thus, we extracted from the simulation a cluster 
with comparable size. 



7.2 Viscosity measurement 

Comparison with experimental viscosity measurement is not an easy task. Viscosity 
measurements are usually performed after solidification under shear rate J8]. Both the 
solidification and the aggregation processes takes place simultaneously. 

In our approach a suspension of dispersed individual globular particles is submitted 
to shear until it reaches its steady state of aggregation. This procedure decorrelates 
solidification and aggregation. However, the solidification process is assumed to have a 
slighter effect than the shear on both the structure and the viscosity. We then compare the 
simulation results with experimental viscosity measurements performed on Al-6.5wt%Si 
alloys partially solidified under stirring ||. Those experimental results were obtained 
with Couette rheometer and the measurements were performed at steady state: when 
the shear stress is stable over time. This experimental procedure allows us to consider 
that the aggregation degree of the microstructure corresponds to the steady aggregation 
degree. The good agreement between simulation results and experimental results plotted 



in figure 16 tends to validate both the former assumption and the basis of the simulation. 



7.3 Particle size distribution 

In figure [17|, we plotted the simulated distribution master curve (see section ||) and ex- 
perimental normalized size distributions resulting from microstructure analysis Jl9| . The 
master curve seem to fit accurately the real distribution profile. 



8 Conclusion 

The structure of semi-solid slurries results from the competition between aggregation due 
to surface forces and break-up induced by the shear field. Aggregation probability is 
set to be constant. Break-up probability is then set to be proportional to the ratio of 
the available kinetic energy over the rupture energy. Using a one-cluster algorithm, we 



supposed that the time evolution of a single cluster approximates the ensemble distribution 
of a population of clusters ("ergodic" hypothesis). Simulation of the time evolution of our 
representative cluster describes accurately the structure of the slurry with no adjustable 
parameter and estimates successfully the viscosity as a function of the applied shear rate 
and the volumic solid fraction. The present work allows a better understanding of the 
structure evolution of sheared semi-solid metallic alloys in the ($, 7) plane. In a future 
work, ripening could be also taken into account through a reinforcement of the necks 
along time after the particles encounter. 

The present approach could be easily transferable to any type of suspension as far as 
the interacting forces between solid particles can be estimated. 

The evaluation of the entrapped liquid fraction could be performed using a general 
dimension code for computing convex hulls. "QuickHulls" from the Geometry Center, 
University of Minnesota, USA, could provide us with a powerful tool. 
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Symbol Meaning 



A Constant depending on thermophysical parameters (1.7 x 10 16 s m~ 3 ) [11 

a Initial particles Radius 

D Diffusion coefficient of the solute in the liquid (5 x 10~ 9 m 2 s _1 ) 

di,d 2 , g?3 Distance between a liquid voxel and the nearest solid particle in the direc- 
tions (Ox), (Oy), (Oz) 

e, e Number and mean number of entrapped liquid voxels in a cluster 

E b Rupture energy of a cluster 

E c Kinetic energy of a cluster 

f c Collision frequency of clusters in a shear field 



K Physical constant defined in equation [10 

M, N 3D Matrice used to store the cluster 

n Number of necks in the fracture plane 

P r Breaking probability 

P a Aggregation probability 

q e Probability for two particles to encounter with favorable crystallographic 
orientation 

R, R Cluster gyration radius and cluster mean gyration radius 

s, s Number and mean number of solid particles in a cluster 

V Volume of the cluster (solid + entrapped liquid) 

x Neck radius 

AL Size of the cluster in the normal direction of the rupture plane 

Act Surface Energy due to liquid / solid surface creation (0.34 J m~ 2 ) 

r] Viscosity of the suspension 

T] Viscosity of the liquid (20 mPa.s) [Q 

7 Shear rate 

7 g Gelation shear rate 

7 S Shear rate from which no aggregation occurs 

$ Volumic solid fraction 

$ e // Effective Volume Fraction (EVF) 

$m Maximum packing fraction 

p Density (2350 kg m" 3 ) 

osl Liquid / Solid surface energy (0.17 J m~ 2 ) 

o h Grain boundary surface energy (negligible: Coincidence Site Lattice hy- 
pothesis) 

r c Contact time between two particles in a shear field 

Tf Formation time of a neck between two particles 

Table 1: Symbols and values for Al-6.5wt%Si 




Figure 1: Rheofluidization: Viscosity decreases between the shear rates % and 1- 
Percolating network, 2- dispersed suspension, 3- suspension of clusters. 




Figure 2: Structure of the suspension in the (7, $) phase-space. 1-Gel, 2-Suspension of 
individual particles, 3-Suspension of clusters, 4-Compact arrangement. In this paper we 
are mostly interested in regions 2 and 3. 




Figure 3: Neck growth between two solid particles. 
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Aggregation 
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Rupture 



Figure 4: Aggregation process: (Al) Initial cluster - (A2) Initial cluster is duplicated - 
(A3) Random rotation - (A4) Sticking of the two clusters. Rupture process: (Rl) Initial 
cluster - (R2) Random rotation - (R3) Random selection of a fracture plane - (R4) Erasing 
the particles in the fracture plane - (R5) Random selection of the resulting cluster. Note 
that the fracture probability depends on the number of neck to be broken in the fracture 
plane. 



Randomly rotate M 



Calculate the break-up and aggregation 
probabilities: P r and P a 
and select whether breaking or sticking occurs 




Delete the solid particles 
in the breaking plane 



Copy M\o N 
and rotate M and N 



Select randomly 
one of the resulting 
clusters, and put it in M 



Stick /Wand N and 
put the result in M 



For each liquid voxel of M, 
select whether it is entrapped or not 



Calculate the new: 

• mean gyration radius: R 

• mean solid particle number: s 

• mean entrapped liquid particle number: e 



Figure 5: Simulation algorithm: the cluster is stored in M and N: 3D matrices. 




Figure 6: Representation of the entrapped liquid (grey) for a 2D Cluster: The liquid is 
entrapped if three or four of the principal directions point to any solid particle (black). 
In the case of two directions pointing to solid particles, the mean distance to solid is 
compared with the cluster gyration radius. 
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Figure 7: Time evolution of the cluster gyration radius. R fluctuates a lot, but its time 
evolution is assumed to be representative of the spatial distribution of the clusters. 




Figure 8: distribution of R after various computation steps. Steady state is reached when 
distribution fluctuations are less than 10%. 



Y= 500 S" 1 



Y= 1 oo s- 1 



Figure 9: Clusters for two different shear rates. Liquid is transparent, solid is grey, and 
entrapped liquid is black. Larger cluster corresponds clearly to higher entrapped liquid 
fraction. 




Figure 10: Gyration radius distribution for 5000 simulation steps at various shear rates 
compared with the mathematical model of Takajo [^3] 




Figure 11: The 7 dependence of the mean gyration radius. For 7 higher than j s , the 
shear breaks the clusters into individual particles. 




Figure 12: The ratio $ e ///$ as a function of 7 for two values of the simulation constant 
corresponding to two initial particle sizes. For 7 higher then the structure is too small 
to trap any liquid. 




Figure 13: Suspension structure depending on 7 and $. 7„ is now $ dependent, except 
for low solid fraction where gelation occurs when a cluster reach the size of the apparatus 
gap (left vertical line). 
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Figure 14: Viscosity as a function of the solid fraction for different shear rate. Higher is 
the shear rate, higher is the solid fraction leading to gellation. 




Figure 15: Experimental and simulated cluster structure taken from an Al-6.5wt%Si semi- 
solid slurry sheared at 900 s~ 1 . (a) 3D reconstruction extracted from reference || (b) 
Simulated with K = 5 x 10 s 
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Figure 16: Comparison between simulation results and experimental viscosity measure- 
ment of an Al — 6.5wt%Si alloy || as a function of the shear rate and solid fraction. 
Simulation results are spread over a domain limited by two reasonable values of a$. This 
domain contains the experimental steady state measurements. 





Figure 17: Normalized cluster size distribution: comparison between experimental mi- 
crostructure analysis of Wan 



19| , and the present simulation. Result of Takajo |24| are 



plotted for comparison. The simulated distribution seem to fit accurately the real distri- 
bution profile. 



